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ABSTRACT 

The operation of an interface dynamo (as has been suggested for the Sun and other 
stars with convective envelopes) relies crucially upon the effective transport of mag- 
netic flux between two spatially disjoint generation regions. In the simplest models 
communication between the two regions is achieved solely by diffusion. Here we in- 
O ' corporate a highly simplified anisotropic transport mechanism in order to model the 

^— > , net effect of flux conveyance by magnetic pumping and by magnetic buoyancy. We 

^ ' investigate the influence of this mechanism on the efhciency of kinematic dynamo ac- 

tion. It is found that the effect of flux transport on the efhciency of the dynamo is 
dependent upon the spatial profile of the transport. Typically, transport hinders the 
onset of dynamo action and increases the frequency of the dynamo waves. However, in 
' certain cases there exists a preferred magnitude of transport for which dynamo action 

Q^ , is most efficient. Furthermore, we demonstrate the importance of the imposition of 

boundary conditions in drawing conclusions on the role of transport. 
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1 INTRODUCTION 



, The solar magnetic field is maintained by the action of a hydromagnetic dynamo. The well-known eleven year solar activity 
5^ ' cycle corresponds to a twenty-two year magnetic cycle, the origin of which lies deep within the solar interior (see, for example, 
5^ , Charbonneau 2005). The key to understanding the generation of the solar magnetic field lies in establishing the form of both 
the large- and small-scale fiuid fiows in the interior. Helioseismology, the study of acoustic fluctuations within the Sun, has 
provided a map of the solar differential rotation (Schou et al. 1998), which shows that the surface latitudinal dependence 
is largely maintained throughout the convection zone while the solar interior rotates as a solid body. These two profiles are 
matched in a thin, stably stratified layer of strong differential rotation (shear) at the base of the convection zone - a region 
known as the tachocline. A discussion of the dynamics of the tachocline can be found in Tobias (2005) and Hughes, Rosner 
& Weiss (2007). 

This region of pronounced shear is believed to be of paramount importance in the generation (via the so-called cu-eSect), 
storage and eventual instability of a strong toroidal magnetic field. The regeneration of poloidal field from toroidal field is 
less well understood. It is typically ascribed to the a-effect of mean field electrodynamics, but the precise origin and location 
of this mechanism remain unclear. A number of possibilities have been suggested (see the review of Ossendrijver 2003). The 
most natural is to appeal to convective turbulence in the presence of rotation (Parker 1955; Steenbeck, Krause & Radler 
1966). There are however some difficulties for this approach both in the linear (Cattaneo & Hughes 2006, Hughes & Cattaneo 
2008) and nonlinear regimes (see the reviews by Diamond, Hughes & Kim 2005; Brandenburg & Subramanian 2005). An 
alternative is to appeal to magnetic buoyancy instabilities of the tachocline field interacting with rotation (Schmitt, Schiissler 
& Ferriz-Mas 1996; Thelen 2000; see also Cline, Brummell & Cattaneo 2003 for a related model). A third possibility places 
the generation of poloidal field well away from the tachocline, the assumed location for the generation of toroidal field. In this 
paradigm (Dikpati & Charbonneau 1999; Dikpati & Oilman 2001) poloidal fiux is regenerated via the decay of active regions 
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Figure 1. Geometry of the two models; (a) Model I, (b) Model II. We consider a local Cartesian geometry in the northern hemisphere, 
vsrith the z-axis pointing radially outwards (upwards), the x-axis pointing towards the equator and the y-axis in the aaimuthal direction. 
The models diflfer only by the 2-extent of the domain. 



at the solar surface (Babcock 1961; Leighton 1969). In these models meridional circulation is invoked in order to transport 
the new poloidal flux to the tachocline. This model too is not without difficulties (sec the review by Tobias & Weiss 2007). 

In this paper we shall consider models of dynamos operating at an interface between regions of turbulent rotating 
convection and velocity shear. The idea of an interface dynamo — with neighbouring regions of a- and w-effects — was 
introduced by Parker (1993) expressly to avoid the difficulties caused by the strong suppression of turbulent magnetic diffusion 
by a weak mean field. In Parker's model, strong toroidal field is generated in the tachocline by radial shear. There, and only 
there, does it suppress the difTusion, leading naturally to confinement of the toroidal field within that layer. The a-effect and 
a healthy turbulent diffusion operate relatively unimpeded in the bulk of the convection zone. Parker suggested that provided 
the diffusivity is not suppressed too strongly, observed solar fluxes and the 22-year solar period could be accommodated. 
Parker's original kinematic model has been extended into the nonlinear regime to investigate explicitly the quenching not 
only of the turbulent diffusion but also of the a-effect (Tobias 1996; Charbonneau &c MacGregor 1997; Markiel & Thomas 



A crucial ingredient for efficient interface dynamo action is the transport of magnetic flux between the two spatially 
disjoint generation regions. For this reason, dynamo action is, typically, more difficult to excite in such models than in those 
with spatially-coincident generation mechanisms (Parker 1993). Indeed, if such a transport mechanism is absent, this places a 
severe constraint on the applicability of interface dynamos (Dikpati et al. 2005). Within the Sun there are however a number 
of possible mechanisms that could transport magnetic field. For example, convective pumping (e.g. Drobyshcvski & Yuferev 
1974) is believed to transport fiux downwards from the convection zone into the tachocline (e.g. Tobias et al. 1998, 2001; Dorch 
& Nordlund 2001). Turbulent diffusion, which is believed to act throughout the convection zone but could be significantly 
reduced in the tachocline (Parker 1993), acts to reduce gradients in field. Mean fiows, such as the meridional circulation 
(Choudhuri et al. 1995), which is observed to be polewards at and just below the solar surface but is undetermined at the base 
of the convection zone, will act to advect mean flux. The only transport mechanism that is unambiguously directed outwards 
is magnetic buoyancy, which is believed to play a key role in the emergence of active regions at the surface (see, for example, 
the reviews by Schiissler 2005 and Hughes 2007). These competing effects may occur on different timescales and depend on 
the local strength of magnetic field. Indeed, the general picture we have of the solar dynamo is of the field being confined and 
maintained at the base of the convection zone, with occasional eruptions. 

In this paper we consider the effect of incorporating non-diffusive flux transport in an interface dynamo model. In keeping 
with the simple yet illustrative nature of Parker's original model, we include magnetic transport via an effective velocity. 
Indeed, such a velocity arises naturally in mean field electrodynamics as the '7-effect', resulting from the anti-symmetric part 
of the a-tensor (Radler 1968). It should be noted that mean field electrodynamics is not without problems, particularly in 
the high magnetic Reynolds number regime applicable to the Sun. Nonetheless, our hope is that it is possible to capture some 
features of the interaction of generation, transport and diffusion of magnetic fields via a simple mean-field model. 

In this paper we shall investigate the role of the 7-effect on the onset of the dynamo instability for two kinematic aoj- 
dynamo models. A natural starting point is to extend Parker's 1993 model, in which the interface separates two spatially 
unbounded regions. It is though important to consider the more realistic case of vertically bounded domains, which is the aim 
of our second model. In §2 we discuss those characteristics that are common to both models. The effect of introducing fiux 
transport is investigated, for the two models, in §3 and §4. The possible implications for an interface dynamo in the Sun are 
discussed in §5. 
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2 FORMULATION OF THE PARKER INTERFACE DYNAMO 

We consider the role of flux transport for two models of kinematic, atj-interface dynamos in a local Cartesian geometry (see 
Figure The region z > represents the convection zone (although there is no convection as such in our framework), wherein 
the a-effect generates the poloidal field and there is no oj-effect. The region z < models the tachocline, where the radial 
shear regenerates the toroidal field (the w-effect) and a — 0. The assumption of axisymmetry translates to solutions that are 
independent of the azimuthal direction, y say, and we shall consider fields that are periodic in x. 
The induction equation for the mean field B takes the form (see, for example, Moffatt 1978) 

^=Vx (C7xB)+Vx£: + 7;oV'B, (1) 

where U represents the mean flow, rjo is the molecular magnetic diffusivity, and £ is the mean electromotive force, which 
takes the form 

= a,,Bj+Aifc|^ + .... (2) 

For simplicity we shall consider the efi'ects of anisotropy only in the aij tensor, which we assume to take the form atj = 
aSij — eijk^k- We assume that Pijk ~ P^ijk', P is then a turbulent magnetic diffusivity. Equation ((!]) now becomes 

^=Vx([7xB)+Vx(7xB)+Vx [aB) - V x {rj {V x B)) , (3) 

in which rj — f3 + rjo. We shall consider models in which rj may take different constant values in z > and z < 0; thus, for the 
purpose of deriving the jump relations below, we do not take rj outside the derivative in the above equation. 

The large-scale velocity U is chosen to represent the radial shear of the tachocline. The effective velocity 7 describes the 
radial non-diffusive transport of magnetic flux between the two layers. We therefore take 

U = U{z)y, 7 = 7^^, (4) 

where dU/dz — GH{—z), 7(2) — j2 + (71 — "/2)H{z) and H{z) is the Heaviside function. Here G, 71 and 72 are constants. 
We may think of the situation in which 7^ < (i = 1, 2) as the case in which magnetic pumping is the prevailing transport 
mechanism leading to a net transport radially inwards, while the case in which the dominant transport effect is by magnetic 
buoyancy (and hence directed radially outwards) is modelled by 7^ > 0. 

As the system is axisymmetric (t/- independent), it is useful to decompose the fleld into poloidal and toroidal parts, 
B{x,z,t) = V X Ay + By. We proceed with the convention that subscripts 1 and 2 refer to quantities in z > and 
z < 0, respectively. Assuming that the shear is much stronger than the a-effect in generating the toroidal field (the acu- 
approximation) , it follows from equation (|3]) that in 2: > the governing equations are 

while in 2: < we have 

dA2 ^2, dA2 dB2 ^2r ^^9^2 dB2 

_^^,V^2-72^, -^-^^^ B2 + G—-J2^. (6) 
In both regions we seek wave-like solutions of the form 

^ = 3fJ[a(z)e*'*+"="], B = K[6(2)e*'*+"="], (7) 

where p = a + icj {a, lu £ k is the wavenumber (fc > 0), and 0(2) and b{z) are the (complex) amplitudes of the dynamo 
waves. Temporally growing solutions (cr > 0), with period T = 2tt/uj, propagate equatorwards when u < and towards the 
north pole when to > 0. Both A and B are continuous. The solutions are matched at the interface 2 = by the jump conditions 
(derived by integrating equation across the interface) 

where [--Jo denotes the jump in the specified quantity across z — 0. Applying the boundary conditions in z, which will be 
described later, leads to a dispersion relation connecting the growth rate and wavenumber to the non-dimensional dynamo 
number {D, say), which measures the relative strength of induction to diffusion. As the definition of D differs between the 
two models, we postpone its definition to the later sections. 

To set the scene we now describe the spatial form of the transport mechanisms that we implement. As already described, 
the efficiency of an interface dynamo is crucially reliant on flux transport between the two generation layers. In Parker's model 
transport is solely by diffusion, with, in the simplest case, 771 =772. The intention of Parker's model was to show how a strong 
fleld could be confined to a thin layer beneath the convection zone owing to that field suppressing the effective diffusion, 
leading to 772 ^ t/i. We shall therefore investigate the role of the 7-effect in models with uniform diffusion (771 =772) and in 
those with 772/771 «C 1. We will consider the following illustrative cases: 

• Bidirectional transport with uniform diffusion 

We consider first the case in which 71 = —72 = 7 with 7 < 0, modelling pumping in 2 > and magnetic buoyancy in 2 < 0. 
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For simplicity we take r)\ = ??2- Although one might expect stronger transport to increase the efficiency of the dynamo, we 
will show that this is not necessarily the case. We will then compare these results with the situation 7 > 0, which represents 
flux transport away from the interface in both generation layers and hence may naively be expected to yield the least efficient 
dynamo. 

• Unidirectional transport with uniform diffusion 

Next we consider 71 = 72 = 7 and rji = ri2. We may think of 7 < as the typical solar case in which transport inwards by 
convective overshoot is expected to dominate the buoyant rise of magnetic field. 

• Dominant transport in the convection zone 

Here we assume that the net transport is much stronger in z > and that effects due to magnetic buoyancy in z > are small 
in comparison with convective transport. We therefore consider 72 = and investigate the effects of increasing the magnitude 
of 71 < 0. This is motivated by what is believed to be the case at the base of the solar convection zone. We compare the case 
772 = rji with the more relevant case rj2 <^ rji. 

For chosen values of the 7-effect and diffusion coefficients, the dispersion relation is solved to yield a critical dynamo 
number Dc and corresponding critical frequency lJc as functions of the wavenumber. Dynamo action sets in for D > D* = 
miiifc [Dc{k)], with preferred wavenumber k = k* and corresponding frequency ujc{k*). We shall say that transport hinders 
dynamo action if increasing the magnitude of 7 increases the critical dynamo number D* , while if D* is somewhere a 
decreasing function of 7 then transport aids the dynamo instability. The effect of transport on the period of the dynamo wave 
(T = 2n/LJc), reflecting the time taken to reverse the sign of the magnetic field, will also be investigated, as will its role in 
modifying the spatial profile of the magnetic field. 

Finally, it is also of interest to consider the effect of transport on the phase relation between the radial and azimuthal fields. 
Observations indicate that BrB^ < at active latitudes, i.e. the two components have a phase shift </> = tt (see, for example, 
Stix 1976; Yoshimura 1976; Schlichenmaier & Stix 1995). This phase difference has traditionally been difficult to reproduce 
with the simplest atj-dynamo models that rely upon diffusive communication: given the sign of dQ/dr from helioseismology, 
the sign of a required for equatorward migration (a < in the northern hemisphere, a > in the southern) generates 
poloidal field whose sign is inconsistent with the above phase relation (see Parker 1987; Schussler 2005). ft is therefore of 
interest to investigate whether more complex atj-models can reproduce the observed behaviour. We note that a non-zero 
phase difference is inherent in the dynamo equations, even for models with spatially coincident generation regions (see Parker 
1955, for example). For our dynamo models, the phase difference will be found by examination of the temporal profiles of the 
azimuthal and radial fields. For example, at a; = 2 = we find 



where <j) is the phase difference. It is important also to note that the observational evidence for the phase lag arises from 
the behaviour of the two components of field at the solar surface. Since the toroidal field is created in the tachocline at 
approximately O.TRq and is subject to complex processes operating throughout the bulk of the convection zone during its rise 
to the surface, the phase lag is likely to vary with depth. In this investigation we choose to evaluate </!> at the interface of the 
convection zone and tachocline {z = 0). For the model of infinite extent in z (Model I) this is done owing to the lack of any 
other preferred location, whereas in the vertically bounded model (Model II) this is done to minimise the role of boundary 
conditions. 

Before turning to the details of the models, it is interesting to note a common characteristic of both models. As can 
be observed directly from equations ([5]) and (|6]), in the large wavenumber limit the evolution is governed by the diffusive 
term. Thus Dc — > cx3 as A; — > 00. Furthermore, since in this limit the critical dynamo number and corresponding frequency 
are independent of the magnitude of the transport mechanism, all the stability curves for different magnitudes of 71 and 72 
collapse onto that for 71 =72 =0. 

3 MODEL I: SPATIALLY INFINITE DOMAIN 

In this section we solve the dynamo equations ((5)1 and (|6]) in a domain that is unbounded in the 2-direction (see Figure [1^), 
with boundary conditions 



since neither the shear nor the a-effect alone can sustain dynamo action. Following Parker (1993), we seek solutions of the 
form (O with 




(9) 



a{z) — + 0, b(z) as \z\ 00, 



(10) 



61(2) =Cexp[-(5' + iO)z], 



ai{z) = (F + Ez) exp [-{S + iQ)z] 



in z > 0, 



(11) 



and 



02(2) = J exp [{s + iq)z], 



h2{z) = (L + Mz) exp [{s + iq)z\ 



in 2: < 0. 



(12) 
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Here C, S, Q, s, q, a, u>, k are real constants, while we allow E, F, J, L, M to be complex. Since our analysis is linear, there 
exists one arbitrary amplitude, which we take as C. The boundary conditions l|10p require that S > and s > 0. Applying 
continuity and the jump relations (O then leads to the relation (see Appendix A) 



[S + iQ + s + iq] 



{S + iQ) + i/{s + iq)+^^^ 



-72) 



2(S + iQ) + ^ 



2(s + ig)-^ 

V2. 



aG 
rjfk^ 



where 



S + K 



= -TT" ± T^\/7? + 4r?i(a + ic^ + mfc''), 
Irii 2r]i 



1 



s + iq = ± 7-^\/7| + 47/2(0- + iuj + 772^2 



1 



(13) 



(14) 



2j72 2r;2 

and /i^ — rj2/r]i. In order to satisfy the conditions at infinity, the signs in expressions (|14p are to be chosen so that S > and 
s > 0. The term in parentheses on the right-hand side of equation (|13p is a dimensionless measure of the relative strength of 
induction to diffusion and is defined as the dynamo number for this problem, i.e. D = aG/'q\k^ . Without loss of generality, we 
take D > 0. There exist similar solutions for Z? < but they propagate in the opposite direction. Note that, owing to the lack 
of any physical length scale in this problem, the dynamo number is here defined in terms of the wavenumber k (e.g. Parker 
1955). Thus when we consider the effects of transport on the onset of the dynamo instability as a function of wavenumber we 
consider the dependence of 7 on k^Dc = (aG/rj\) ^. 

We note that in the absence of the 7-effect the model reduces to that of Parker (1993). In the case of uniform diffusion, 
(r;i = r]2) equations (|13|l and (|14p then yield 

'aG 



±32fc-' 



±rj\k 



(15) 



In what follows we take the {a.G/'ql)c > 0, ujc > Q solutions as our starting point when considering the effects of non-zero 
transport; i.e., as the magnitude of 7 increases, our solutions are smoothly linked to expressions (|15l) with the positive signs. 
As discussed earlier, we split the general case into the following subcases. 



3.1 Bidirectional transport and uniform diffusion 

We begin with the case 71 = —72 = 7 and ?7i = 772 = r). Relations p3|l and (|14|l then yield 

(^) =^(7^ + 8,V)(7^+4,V)^/^ ^. = |(7^+4r,V)^/^ (16) 

It is clear that increasing 7 increases both {oiG/rji^ ^ and ujc- Thus transport hinders the onset of the dynamo instability and 
decreases the time taken to reverse the field polarity. Furthermore, both relations are independent of the sign of 7. Thus, 
although one may expect the dynamo properties of the two cases 7 < and 7 > 0, representing flux transport towards and 
away from the interface in both domains, respectively, to be very different, the critical parameters governing the onset of 
dynamo action are not affected by the direction of flux transport. We believe this urmsual property to be a consequence of 
the infinite extent of the model. 

It is noted that the parameters governing the z-structure of the solution are sensitive to the sign of the 7-effect, viz. 
S -\- \Q = s -\- iq = —7/27; -|- ■y/72 -I- 4r;(cr -|- ioj -I- r]k'^ ) /2r;. As a consequence, the dependence on 7 of the phase difference 
between the azimuthal and radial fields is different for positive and negative 7, as shown in Figures [5][a,b). In the case when 
transport is directed towards the interface (7 < 0), the phase lag is increased as I7I increases, while transport away from the 
interface (7 > 0) brings the two components of field more in phase. 

The structure of the azimuthal and radial fields is illustrated in Figures [Hc,d). As expected, the eigenfunctions are 
narrower (wider) when transport is towards (away from) the interface. 



3.2 Unidirectional transport and uniform diffusion 

Given that flux transport towards the interface hinders the excitation of the dynamo instability in exactly the same manner as 
does transport away from the interface, it appears intuitively unlikely that a unidirectional 7-effect operating throughout the 
domain will enhance the efficiency of dynamo action. Indeed, in the case 71 = 72 = 7 and 7)1 = 772 = rj, analytic investigations 
yield 

(f)^--.-(i.ji_)', („, 

The dynamo is most easily excited at k = kc — |7|/(27]\/3), where (aG/r;i)^ — 64|7|''/ (Sri^VS). Thus transport again hinders 
the onset of the dynamo instability. The frequency increases with increasing magnitude of the transport, and the dynamo 
prefers to generate waves with ever decreasing wavelength (A = 2n/k). 

It may also be shown analytically that at z = increasing the magnitude or changing the direction of the transport 
has no effect on the phase difference between the toroidal and radial fields, as = 7r/4 always. Furthermore, since the phase 
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Figure 2. Model I: Bidirectional transport and uniform diffusion (71 = —72 = 7, f?2 = '?! = !)■ (a) Phase difference {j>, (b) Phase 
difference divided by the period (/i/T, (c) Toroidal field, (d) Poloidal field. 7 = represents no transport, 7 = —1 transport towards the 
interface and 7 = 1 transport away from the interface. The eigenfunctions are for k = 1 and the phase difference is calculated for A: = 1 
and X = z = 0. 

relation remains constant even though the critical parameter governing the onset of dynamo action can increase dramatically, 
the phase lag is not a property that controls the efficiency of dynamo action in the regime studied here. 

Figure [3] compares the structure of the azimuthal and radial fields for unidirectional transport directed downwards with 
that for upwards transport. Downwards transport reduces the decay length of the fields in z > 0, while upwards transport 
reduces the decay length in z < 0. 

3.3 Dominant transport effect in z > 

We now consider the effects of removing the 7-effect from z < 0. In these cases the toroidal field generated in the shear layer 
will rely upon diffusion to transport it into z > 0, where it will act as a source for poloidal field regeneration. 



3.3.1 Uniform diffusion 

Figures |4ja,b) illustrate the marginal stability curves for the case 71 < 0, 72 = and rji — rj2 = rj = 1. For small wavenumbers 
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^8 
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1/3 
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as fe — » 0, 



(18) 



which may be compared with (I15|l for the case of 7 = 0. Thus, again, transport hinders the onset of dynamo action and 
increases the frequency of the dynamo waves. It is interesting to note that in comparison with the uniform transport case of 
^3.2\ the removal of the 7-effect from z < results in a phase difference that is dependent on the magnitude of 71 [FigurelUc)]. 
In this case increasing 71 increases the phase difference. The ratio of the phase difference to the period of the wave also increases 
[Figure |4jd)]. The eigenfunctions, shown in Figures |4je,f), exhibit similar behaviour to the case of a unidirectional transport, 
with the fields being distorted towards the interface from above (the direction in which the transport operates). 
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Figure 3. Model I: Unidirectional transport and uniform diffusion (71 = 72 = 7, J?! = = li = !)■ (a,b) Toroidal and poloidal fields 
when the 7-effect is directed downwards. (c,d) As for (a,b) with the 7-effect directed upwards. For 7 = —50, 10, 50, the azimuthal field 
has been multiplied by 10, 50, 1000, respectively. 



3.3.2 Reduced diffusion in z < 

The marginal stability curves for the case 71 < 0, 72 = 0, 771 = 1 and 772 = 0.01 are shown in Figures [5ja,b). Again 
transport hinders the operation of the dynamo and decreases the period of the dynamo wave. The phase difference between 
the radial and azimuthal fields behaves in a qualitatively similar manner to the cases investigated above in which the 7-efFect 
is directed downwards, with both (f) and 4)/T increasing as the magnitude of 71 increases [Figures [5l[c,d)]. The behaviour of the 
eigenfunctions appears to be such that the downwards-directed transport reduces the decay lengths in z > [Figures [5je,f)]. 

Comparing these results with the uniform diffusion case described above f i)3.3.1|l . one finds that in the absence of 71 the 
dynamo is easier to excite in the case of uniform diffusion than for ?72/'?i = 0.01, for those wavenumbers considered. Indeed, 
in the limit = ri-z/rji 0, Parker (1993), in the absence of transport, derived the relations 

faG\ ^2^/3 3 ^Vif^'^\2 

WJ. ^ ' '^=^^17^'=' ^^^^ 

while for /i^ = 1 we have expressions (|15p . Clearly, in the absence of non-diffusive transport, /^^ ^ yields a greater critical 
dynamo number and a smaller frequency than for the case of uniform diffusion. However, our results reveal that for finite 
71 the situation is reversed at a certain wavenumber, k — fc(7i), say. For wavenumbers less than k the dynamo is easier 
to excite in the case with smaller 772, while for k > k dynamo action occurs for smaller values of (aG/rji) in the case with 
uniform diffusion. This can be seen by comparing Figures |4j a) and[5{a) for the case 71 = —1, where k{—l) « 0.3. For those 
71 considered, the frequencies are larger in the case of uniform diffusion than they are for 772 = 0.01. However, it should be 
noted that for different spatial profiles of 7 it is not always the case that uniform diffusion results in a larger frequency (for 
example, compare the frequencies in Figure IUb) for 772 = 1 and 772 = 0.01 in the case of 71 =72 = —1)- Finally, inspection 
of the eigenfunctions shows that the reduced 772 results in a smaller vertical length scale for the toroidal field in z < 0, as 
expected (Parker 1993; Tobias 1996). 
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Figure 4. Model I: Dominant transport in z > and uniform diffusion (71 < 0, 72 = 0, ri2 = Vi = !)• (^) Critical parameter k^Dc, (b) 
Critical frequency lo^, (a) Phase difference <j>, (b) Phase difference divided by the period <j>/T, (c) Toroidal field, (d) Poloidal field. The 
eigenfunctions are for k = 1 and the phase difference is calculated for fc = 1 and x = z = 0. 



3.4 The competition between the 7-effect and diffusion 

It is instructive to consider the effect of varying the diffusivity in 2; < in the absence of the 7-effect and to compare the 
results with the corresponding situation when 7 is non-zero. We consider a fixed value of the wavenumber, here chosen to be 
fc = 1, and we take 71 = 72 = 7 for simplicity. Figure |6j a) illustrates that in both cases the effect on {aG/rjf)c of varying rj2 
is non-monotonic. Decreasing r]2 from unity initially allows a more effective production of toroidal field in z < and thus the 
dynamo becomes easier to excite. However, there exists a transition point {ri2{-f), say) after which decreasing the diffusion 
hampers the coupling between the two regions sufficiently that the onset of dynamo action is delayed. 

Figure [6jb) illustrates that when the 7-effect is absent the frequency decreases monotonically as rj2 decreases, while when 
7 7^ there are two distinct regimes. The first is when 7 <^ kr]2, in which the frequency of the wave is governed by diffusive 
transport between the layers. Here the frequency decreases as 772 decreases. The second regime begins when kr]2 becomes 
comparable with 7, and the frequency is then controlled both by diffusion and the 7-effect. As 772 is decreased further there 
is an increase in the frequency of the dynamo wave until ujc eventually saturates at a constant value. In the final state the 
frequency of the wave is governed by the 7-effect. 

Figure [6jc) shows that the phase difference between the toroidal and radial fields is also governed both by diffusion and 
the 7-effect. The phase difference is typically a non-monotonic function of the diffusion, even for the case 7 = 0. The two 
components of field initially become more out of phase as the reduced diffusivity favours the generation of toroidal field in 
z < 0. However, as the diffusion is reduced further the phase difference decreases. The ratio of the phase difference to the 
period of the wave is shown in Figure [6jd). 

The effect of reducing the diffusivity in 2: < on the spatial structure of the fields is illustrated in Figures [6je,f). In the 
case 7 = (compare Figure |3| a) with[5je)) the effect of decreasing r]2 is increasingly to confine the toroidal field beneath 
2; = 0, as first shown by Parker (1993), while when 7 7^ a sufficiently small diffusion in z < allows significant amounts of 
field to be advected away from the interface. 
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4 MODEL II: SPATIALLY CONFINED DOMAIN 

In this section we investigate a more realistic model whose vertical extent is limited to 2 £ [— L,L], where L is finite (see 
Figure [T]d). All other aspects of the geometry remain as in Model I. One consequence of the finite domain is that we may 
meaningfully non-dimensionalise the system. We let 



[x,a, G = ^,r;,7,B,yl,t] = [Lx', aoa', GoC, r;o?)', 707', -BoS', — f'], 

where primed variables are dimensionless. The dynamo equations then become (dropping the primes) 

—— ^ aB + riiV Ai -Yi——, -7— ^ -qiV Bi -Yi—— m 2 > 0, 
at oz at az 

and 

in z <Q. 



dA2 ^2 , ^ dA2 



(20) 

(21) 
(22) 



Here Yi,2 = R-t'ii,2 is a dimensionless measure of the 7-effect and R-y = 'yoL/rjo is its Reynolds number. The dynamo number 
is defined as D = aoGoL^ /rjo^ > 0. 

Travelling wave solutions of the form ([7} are sought in each of 2: > and z < and the solutions are matched at the 
interface z — Ohy the jump relations ((8}. The boundary conditions that approximate the matching onto an external potential 
field (see Zeldovich, Ruzmaikin & SokolofT 1983) are, 



Bi{x,z^l,t)^0, —l(x,z = l,t)^0, 



B2{x,z = -l,t) = 0, 



dA2, 
-^(x,z ■ 
oz 



^l,t) =0. 



(23) 
(24) 



The resulting dispersion relation, along with details of its derivation, can be found in Appendix B. Again we split the general 
case into a number of subcases. 
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Figure 6. Model I: The competition between the 7-effect and diffusion (71 = 72 = 7 < 0, »7i = 1, fc = 1). (a) Critical parameter k^Dc, 
(b) Critical frequency lUc, (c) Phase difference 0, (d) Phase difference divided by the period (p/T, (e) Toroidal field, (f) Poloidal field. 
The eigenfunctions are for 7 = —0.05; for t]2 = 0.01 and 772 = 0.001 the toroidal field has been divided by 10 and 20, respectively. The 
phase difference is calculated for x = z = 0. 



4.1 Bidirectional transport and uniform diffusion 

It is instructive first to investigate the case in which the 7-efrect is directed towards the interface from both above and below. 
We therefore take Yi — —¥2 — Y with Y < and for simplicity we consider the case of uniform diffusion, J71 = J72 = V- 
The marginal stability curves and corresponding frequency profiles for different magnitudes of the transport are shown in 
Figures [7ja,b). For fc ^ 1, all profiles collapse onto 

Dc«32(^^^, LOc^r,k\ (25) 

as the system reduces to Model I in the large wavenumber limit (cf. (jlSp : Parker 1993). The numerical results illustrate that 
the dynamo instability is easiest to excite ai k <^ 1, where ~ Do{Y,ri) and ujc ~ uji{Y,rj)k. Analytic investigations for 
y = yield Do(0, 77) — IBjyoJi and uJi{0, rj) — r\\fzfl. Figures [7j a, b) illustrate that increasing the magnitude of the transport 
decreases both Do and the corresponding frequency of the dynamo wave. Thus here the transport effect aids the onset of the 
dynamo instability. The time taken for the field to reverse sign increases. 

It is interesting that both of the above results concerning the dynamo number and frequency are the direct opposite 
of any situation that we have previously encountered, i.e. the critical dynamo number and frequency in ^ were increasing 
functions of 7. In particular, this is true for the corresponding case of Model I (see i^3.1|l . and we are led to the conclusion 
that boundaries in the z-direction play a crucial role when considering the effects of transport on the efficiency of dynamo 
action. In fact, investigating the complementary situation in which 7 is directed away from the interface in both domains 
(Vi = —Y2 = Y,Y>0) yields numerical results with coefficients Do and loi monotonically increasing functions of Y. Recalling 
that Model I was dependent only on the magnitude of the bidirectional transport mechanism, we conclude that the boundaries 
in the z-direction cause the model to become sensitive to the direction of flux transport, which is more realistic. 

Figures[7jc,d) illustrate the effect of transport on the phase difference between the azimuthal and radial fields, and on the 
ratio of the phase difference and the period of the wave. For comparison, also presented is the case when transport acts away 
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Figure 7. Model II: Bidirectional transport with uniform diffusion (Yi = —Y2 = Y, K < 0, r;i = 772 = 1). (a) Critical dynamo number 
Dc, (b) Critical frequency tJc, (c) Phase difference </), (d) Phase difference divided by the period fli/T, (e) Toroidal field, (f) Poloidal field. 
The toroidal field has been divided by factors of 2 and 4 for Y = —5 and Y = —10 respectively (fc = 1). The phase difference and 0/r 
are compared with the case of y > representing transport away from the interface in both z > and z < (A; = 1, a; = z = 0). 



from the interface in both regions. Transporting flux towards the interface causes the two components of field to become more 
out of phase whereas when Y is directed away from the interface the phase difference decreases, in agreement with Model I. 

Conveying magnetic flux towards the interface favours fleld generation there, causing both the toroidal and poloidal 
fields to become increasingly localised around 2; = [Figures [7je,f )] . In the contrasting situation of flux transport away from 
the interface {Y > 0) the main effect is to cause the peak in toroidal magnetic field to move towards the lower boundary 
(behaviour that is qualitatively similar to that shown in Figure ISje) for uniform downwards transport), whereas there is little 
change in the z-structure of the poloidal field. 



4.2 Unidirectional transport and uniform diffusion 

When the transport velocity is uniform throughout the whole domain {Y\ = Y-i = Y^ and the diffusivities are equal, the 
dispersion relation is independent of the sign of Y . Figures [Ha,b) illustrate that for k <^ \ the critical parameters exhibit 
different behaviours depending upon the magnitude of the transport: for 1^1 < Ya(r]), say, the most readily destabilised mode 
satisfies Dc ~ i3o(y, ??), i^c ~ ijo\(Y^r\)k, whereas for |y| > Ys we have Dc ~ D_i(y, 7))fc~^, ojc ~ {joq^Y^T]) [for 77 = 1, Ys ~ 3.1 ]. 
Since at large wavenumbers the marginal stability profiles must approach the results (|25[) for Model I, it follows that for those 
|y| > Ys there exists a preferred wavelength for the dynamo instability (for example. A* ~ 27r/3 for Y — —10). Figures [8ja,b) 
illustrate that transport hinders the excitation of the dynamo and decreases the period of the wave. 

The phase difference between the azimuthal and radial fields follows a non-monotonic behaviour as the magnitude of Y 
is increased, and is dependent upon the direction of flux transport [Figures [8|c,d)]. When y < the two fields are most in 
phase for Y « —1.6, which can be compared with the situation for Y > where the stationary point is a maximum and the 
fields are most out of phase for Y ~ 3.2. 

The spatial structure of the toroidal and poloidal fields in the presence of a uniform transport velocity directed towards 
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Figure 8. As for Figure [7] with Yi = Y2 = Y < 0, rji = 1^2 = 1- The toroidal field has been divided by factors of 10 and 50 for Y = —5 
and Y = — 10 respectively. 



the boundary z — —L is shown in Figures IHe,f). Interestingly, increasing the magnitude of the transport results in the 
majority of the poloidal field residing in ^ < 0, even though it is generated in 2; > 0. 



4.3 Dominant transport effect in 2 > 

We now consider the case motivated by the Sun, in which the downwards pumping of the convection is much stronger in z > 
than in 2 < 0. We take Y2 = and Yi < and consider the effects of increasing |Yi|. Flux transport in the tachocline region 
is then controlled only by diffusion. 



4-3.1 Uniform diffusion 

We first consider the case of uniform diffusion, ?7i = 772 = Figures |9ja,b) show that there is a magnitude of the transport 
(¥1^(77), say) for which there is a qualitative change in the behaviour of the solution at small wavenumbers. This behaviour 
is similar to the case of a unidirectional 7-effect (see § I4.2|l although here the results are dependent on the direction of fiux 
transport. For rj — 1 and Yi < we have Yi^{l) « —4.9. For |Yi| < the most readily destabilised mode satisfies 

Dc ~ Do{Yi,r]) and ljc ~ uji{Yi,ri)k for fc <C 1, while for |Yi| > and fc ^ 1 we have Dc ~ D-i{Yi,ri)k~^ and 

u!c ~ uJo{Yi,rj). Most notably however, and in contrast to all previous calculations, it is found that the minimum of Dc over 
all wavenumbers as a function of the transport speed Yi possesses a stationary value. Thus there exists a preferred magnitude 
of transport for which dynamo action is most efficient. For rj = 1 the transport that leads to the most efficient dynamo is 
Yi* ~ —1.2, where Dc ~ 16.8. We shall return to this point in MA\ [see also Figure [TU cI] . 

The effect of the transport mechanism on the phase relation between the radial and azimuthal fields is illustrated in 
Figures |9jc,d). The monotonically increasing behaviour with |Yi| is similar to that for the same case in Model I (cf. Fig- 
ures |4jc,d)). 

Figures IUe,f) illustrate the structure of the poloidal and toroidal fields as the transition from transport enhancing to 
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Figure 9. As for Figure [7] with Yi < 0, Y2 = 0, Vi = V2 = 1- The toroidal field has been divided by factors of 2 and 4 for Y = —5 and 
y = — 10 respectively. 



hindering the dynamo occurs. Although it is difficult to infer from the eigenfunctions alone why a small magnitude of transport 
helps the dynamo (there appears to be little difference in the spatial profile of the solutions in the cases Yi =0,-1), it is 
interesting to note that in the case when the transport exceeds the preferred magnitude the majority of both the toroidal and 
poloidal fields resides in 2 < 0. This behaviour is suggestive that corresponds to the strength below which the main effect is 
to enhance the downwards transport of flux into the tachocline, more readily feeding the u-eSect and leading to more efficient 
dynamo action, whereas above ¥{ the diffusive transport of toroidal field into the convection zone is hindered significantly, 
thus inhibiting the operation of the a-effect and making the dynamo cycle more difficult to complete. Interestingly, the 
behaviour of the fields in z > bears some similarity to the case of a uniform downwards transport throughout the domain 
[cf. Figures |8jc,d)]. Thus the reason that the dynamo is more difficult to excite in the case of uniform transport (the minima 
of Dc over all k in Figure |8ja) are greater than those in Figure |9l a)) is likely to be due to the behaviour in 2 < 0: when Y2 is 
finite the toroidal field is translated towards the lower boundary where it is annihilated. 



4-3.2 Reduced diffusion in z < 

We now consider the case in which the the diffusivity in 2 < is much smaller than that in 2 > 0. We take 771 — 1, i]2 ~ 0.01 
and compare the results with the uniform diffusion case above f i]4.3.1|l . The profiles of the critical dynamo number and 
corresponding frequency are shown in Figures [TOT a.b). 

At large values of the wavenumber the profiles collapse onto those for Vi = I2 = 0. Indeed, for k ^ I and /i^ ^ 
the critical dynamo number and frequency are given by the relations (|19|) (Parker 1993). The numerical results show that 
for A: ^ 1, Dc ~ Do{Yi,rii,ri2) and luc ~ aji(Yi, 772)^, where the coefficients Do and uji are non-monotonic in Yi. Thus 
again a preferred magnitude of transport, ¥{ , exists for which the dynamo is most easily excited. For 772 ~ 0.01 and Yi < 0, 
Yi ~ —6.5 [see also Figure [TTT d)]. It is noted that when transport is absent the dynamo is easier to excite when 772 = 1 than 
when 772 = 0.01; however a comparison of the preferred strengths in the two cases (i.e. Dc{Yi ~ —6.5,772 = 0.01) — 0.99 
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Figure 10. As for Figure [7] with Yi <Q,Y2= 0, 771 = 1, rj2 = 0.01. 



compared with DdY^ ~ —1-2, J?2 = 1) = 16.8) shows that the stronger 7-effect and smaller diffusivity ratio results in a more 
efficient dynamo. 

The effect of transport on the phase relation between the radial and azimuthal fields, and also the ratio of the phase lag 
to the period of the wave, is shown in Figures [TOT cd). The behaviour is qualitatively similar to that for the case of uniform 
diffusion with Y2 — Q (cf. Figures |9lc,d)). 

Figures lior e.f) illustrate the effect on the eigenfunctions of increasing the transport in 2 > 0. The structure of the 
toroidal field is affected only weakly, whereas the poloidal field becomes more localised around the interface, the effect being 
particularly pronounced in 2: > 0. 



4.4 The most efficient dynamo 

Finally, to summarise the findings for this model and to compare the efficiency of the four cases considered. Figures ITIT a-d') 
present the minimum value of the critical dynamo number over all wavenumbers as a function of the 7-effect for each case. 
Figure lllf a) illustrates the most efficient dynamo model in which the 7-effect is directed towards the interface from both 
above and below. Increasing the strength of the 7-effect in this case monotonically decreases the threshold for the dynamo 
instability. Figure lllf b') illustrates the least efficient dynamo model. In this case a uniform transport mechanism operating 
throughout the domain hinders dynamo action, with the minimum critical dynamo number increasing monotonically with 
increasing strength of the 7-effect. Figures [TlT c.d'l illustrate the cases when the non-diffusive transport is restricted to operate 
only in 2 > 0, with uniform diffusion in the former and a reduced diffusion in z < in the latter. In both of these cases there 
exists a preferred magnitude of the 7-effect for which dynamo action is most efficient. 
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Figure 11. Model II: minfe[Dc(fc)] versus magnitude of transport, (a) Bidirectional transport towards the interface with uniform diffusion 
(Yi = —Y2 = Y , Y < 0, rji = rj2 = 1), (b) Unidirectional transport and uniform diffusion (Yi = Y2 = Y , rji = ri2 = 1), (c) Dominant 
transport in 2 > with uniform diffusion {Yi < 0, Y2 = 0, '?i = »?2 = 1), (d) Dominant transport in 2 > with reduced diffusion in 2 < 
(Yi < 0, 1^2 = 0, rji = 1, J72 = 0.01). Cases (c) and (d) have a finite most efficient strength of the 7-effect, Yj* —1-2, —6.5, respectively. 



5 CONCLUSION 

We have considered the role of anisotropic transport on the efficiency of an interface dynamo in the linear regime, and its 
effect in modifying the temporal structure of the solution, through the frequency of the dynamo wave and the phase relation 
between the radial and azimuthal fields. One of the main findings of the investigation is that the conclusions to be drawn 
depend crucially upon the geometry of the model. 

The first model considered was of infinite vertical extent. In the absence of transport, the preferred wavenumber for 
the onset of the dynamo instability is A; = 0, where the critical parameter {aG/rj\)c ~ 0. Although, clearly, this cannot be 
decreased by incorporating the 7-effect, for all forms of transport investigated it is found that minfc[(aG'/?7i)c] increases as 
the magnitude of the transport increases. Thus dynamo action is less readily achieved. Also, the period of the dynamo wave 
decreases as the magnitude of the 7-efrect increases. Thus although the dynamo is more difficult to excite, it more quickly 
reverses the polarity of the field. Perhaps most noteworthy is that for this model the two cases of the 7-effect directed towards 
the interface from both above and below and of the transport directed away from the interface in both regions have identical 
critical dynamo numbers and corresponding frequencies fi )3.1[) . 

The effect of the transport mechanism in the more realistic Model II is more complex. The model geometry is the same 
as that of Model I except that the domain of the problem is bounded in the vertical direction. Here it is found that a uniform 
transport velocity throughout the convection zone and tachocline delays the onset of the dynamo instability and increases 
the frequency of the dynamo waves. However, transport directed towards the interface in both z > Q and z < Q decreases 
both the critical dynamo number and the frequency. Thus transport aids the onset of the instability in this case, although it 
then takes longer to reverse the field polarity. Restricting the transport mechanism to the convection zone reveals a preferred 
strength for dynamo action. 

It is of interest to relate our studies to previous investigations of the effects of transport on dynamo models. Brandenburg 
et al. (1993) and Gabov et al. (1996) considered models of galactic dynamos and investigated the role of turbulent diamagnetism 
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on the dynamo instability via the additional term V x (— V/9/2 x B) in the induction equation, where /3 is the turbulent 
diffusivity. Brandenburg et al. (1993) illustrated that a galactic wind can enhance the instability by transporting field into 
regions of stronger generation, while turbulent diamagnetism can act to oppose this effect and hence hinder dynamo action. 
Similarly, Gabov et al. (1996) illustrated that the critical dynamo number is larger when the effects of turbulent diamagnetism 
are present. This was attributed, in part, to an increase in the effective diffusion coefficient. It was then shown that the effective 
dynamo number Dc/p (where P is the mean value of the turbulent diffusion coefficient) decreases as /3 is increased. It was 
further illustrated that for highly supercritical dynamo numbers the linear growth rate is larger when turbulent diamagnetism 
is present than in its absence. In our study we have concentrated on the effects of transport on marginally stable solutions. 
Nightingale (1983) incorporated the 7-effect into a spherical a^-dynamo model and found that the dynamo growth rate is 
smaller when 7 is non-zero than when there is no transport. Radler (1986) considered spherical au)- and a^-dynamos and 
showed that although radial transport mainly hinders dynamo action, certain strengths of 7 can promote the instability. 

Finally, we remark that interface dynamo models relying purely upon diffusive communication are sometimes rejected 
owing to their inability to reproduce either the observed direction of propagation of the dynamo belts (exemplified in the 
butterfly diagram) or the observed phase relation between the radial and azimuthal fields. In this regard, it is of interest to note 
that the study by Brandenburg et al. (1992) showed that a sufficiently strong downward pumping can lead to an equatorward 
migration of the magnetic field even when the conventional rule for the sign of adQ./dr yields a poleward migration. Also, 
here we have shown that the transport mechanism changes the ratio of the phase difference to the period of the dynamo wave. 
It is therefore entirely possible that tuning the 7 parameter could enable interface dynamo models to reproduce the observed 
behaviour. 
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APPENDIX A - DISPERSION RELATION FOR MODEL I 

In 2 > we seek solutions of the form 

Bi = Cexp [t{a + iu)] cxp [-z(S + iQ)] oxp {ikx), (26) 

= (F + Ez) exp [t{cr + iuj)] cxp [-z{S + iQ)] exp [ikx), (27) 

where C, E, F are constants and a, uj, S, Q, k are real, with S > so that B — > as z — > +00. Substitution into equations 
(HI gives 

E[27^i{S + iQ)+'yi] = aC, a + icj ^ -k'^j^i + rjiiS + iQf + jiiS + iQ), (28) 

with F arbitrary. 

In z < we look for solutions of the form 

A2 = J exp [t{a + iuj)] exp [z{s + iq)] exp (ikx), (29) 

B2 = {L + Mz) exp [t{(j + icj)] exp [2(5 + i^)] exp {\kx), (30) 

where J, L, M are constants and ct, uj, s, q, k are real with s > so that B ^ as z ^ —00. Substitution into equations ^ 
then gives 

M[2»72(s + ig) -72] = -ifcGJ, a + ioj = -fc^?72 + 772(5 + ig)^ - 72(5 + ig), (31) 

with L arbitrary. Applying continuity of A and B across z = Q, together with the jump conditions yields 

F = J, C = L, 72L - 71C -r;2M- 772^(5 + ig)-r?i (5 + iQ)C = 0, E = F[S + iQ) + J{s + ].q). (32) 

Equations (I28b 'l. (|31b .) and ((32]) may be combined to yield two equations for C and F. A non-trivial solution requires that 
the determinant of the coefficients of C and F vanishes, from which we obtain the relation 

(S + iQ + s + -iq) {{S + iQ) + At'(s + ig) + '-'^'^ '^'-' ^ {2{S + iQ) + ^2(s + iq) - ^) = '^^^D, (33) 

where D = aG/rjik'^ and fi^ = ri2/r)i. Recall that we also have equations p8b ) and (|31b) that relate s, q, S and Q to the 
complex growth rate. 
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APPENDIX B - DISPERSION RELATION FOR MODEL II 

We consider first the convection zone (2 > 0). The dynamo equations (|2H) . together with the proposed form of solutions ([7]), 
yield 

(fbi Yi dbi qf ^ _^ d^ax Yi dai gjai _ aW ^^^^ 
dz^ 771 dz rji dz"^ 771 dz rji 771 

where ql = p + i]ik^ and Fi = R-yji- Solving equation (|34b ,) for bi and applying boundary condition (|23b l yields 

61(2) = 2Asinh[m(2 - 1)] exp ( ^'^f^^ , (35) 



2771 

where ^ is a constant and rTT = \/Y-^ +4771^^/2771. Substituting bi into equation (I34b l and solving the homogeneous equation 



we obtain 



ai^ (2) = exp (^^^^— (Ccosh[7n(2 - 1)] + £sinh[7n(2 - 1)]) . (36) 



Looking for a particular solution of the inhomogeneous part of the form 



ai^(2) = 2exp ( ^^^£ ^ ] (/cosh[m(2 - 1)] + Fsinh[m(2 - 1)]) (37) 



(38) 



2771 

we find I = — QA/77irn and F = 0. On application of boundary condition (|23b) we finally obtain the full solution 

01(2) = exp ( — !-^£ 3^ I ( Ccosh[m(2 — 1)1 + (1 + — i- | sinh[m(2 — 1)1 — cosh [771(2 — 1)1 | 

^' ^\ 27^1 J\ ^ ^ -^^ [7717712^ 2771 ^ 27;imJ ^ m7?i ^ ^ '^J 

Following a similar procedure in the tachocline region (2 < 0), the corresponding equations to (|34p read 

d2&2 Y2 db2 ql , \kDGa2 d^a2 Y2 da2 ql „ 

TTl = , -r^ ^ 02=0, (39) 

a2^ 772 dz 772 772 dz'^ 772 dz 1)2 

where gl = P + 772^^, with solutions 

a2 = Jexp f Y2l^±A\ ( cosh[n(2 + 1)1 - J^sinh[n(2 + 1)1 1 , (40) 

\ 2772 / \ 277277 

■Y2{z+^\ fikDGY2J ur-r^iM^^r '^kPGJ z 

^r^T — 1 + 2)cosh 77 2 + 1) + L 

2772 / V 47;|n2 ^ ' L ^ ^ ;j -r 

where n = ^ Y^ + 4772 72772 . Applying the jump conditions (|8]l (with 7 replaced by Y) and continuity of A and B across 
2 = gives four equations for the four unknown constants A, C, J and L. A non-trivial solution requires that the determinant 
of the coefficients vanishes, from which we obtain the dispersion relation 

—2(27/271 cosh[n] — I2 sinh[n]) j^ioDGVi 1^2^771 cosh^ [m] cosh[n] + 

cosh[77i] sinh[7n] ( - iaDGY^Yik cosh[n] + 477i7n2n(-y2^ + 477271.^) sinh[n]) - 

7nsinh2[m]('[iaDG(277i + yi)y2fc + 4772 (Yi^ - \r]\m)n] cosh[n] 



62(2) = exp ( ' ) ( ^^2-2 (1 + ^)cOsh[77(2 + 1)] + ( L - ) Sinh[77(2 + 1)] ) , (41) 



+2n(77iyiy22 - 7;2V2(yi' - ^-qW) - A!q^r]lY^f^) sinh[7i])] + 

sinh[n] |iaDG(2772 -y2)A:(-2(77 1 + yi)7r7, + 2771777, cosh [2771] + yi sinh[277i]) (— y2 cosh [71] + 277277 sinh[7i]) 
— 4rn,n(27;i7n cosh[7n] + yi sinh[7n]) |^ — 27;|(yi2 — 47/1 rri2)n cosh [n] sinh[7tt] + 
- 277i7n(y2^ - 477277,2) cosh[7Tt] + ( - 77iyiy2^ + r]2Y2{Yl - 4771 7n,2) + A,n^r]lYYn) sinh[7ft]) sinh[n]] | = (42) 
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